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ABSTRACT 



Context. The 6.7 GHz maser transition of methanol has been found exclusively towards massive star forming regions. A majority 

of the masers have been found to lack the presence of any associated radio continuum. This could be due to the maser emission 

originating prior to the formation of an Hii region around the central star, or from the central object being too cool to produce a Hii 

region. 

Aims. One way to distinguish between the two scenarios is to determine and model the spectral energy distributions (SEDs) of the 

masers. 

Methods. We observed a sample of 20 6.7 GHz methanol masers selected from the blind Arecibo survey, from centimeter to submil- 

limeter wavelengths. We combined our observations with existing data from various Galactic plane surveys to determine SEDs from 

centimeter to near-infrared wavelengths. 

Results. We find that 70% of the masers do not have any associated radio continuum, with the rest of the sources being associated 

with hypercompact and ultracompact H ii regions. Modeling the SEDs shows them to be consistent with rapidly accreting massive 

stars, with accretion rates well above 10~^ M© yr~^ The upper limits on the radio continuum are also consistent with any ionized 

region being confined close to the stellar surface. 

Conclusions. This confirms the paradigm of 6.7 GHz methanol masers being signposts of early phases of massive star formation, 

mostly prior to the formation of a hypercompact H ii region. 

Key words. Masers - Stars: formation - H ii regions 



1. Introduction 

Class II methanol masers, of which the 6.7 GHz maser is the 
strongest, are commonly found towards massive star forming re- 
gions. A number of early surveys for 6.7 GHz methanol masers 
were focused towards IRAS sources wit h colors sug gestive of 
ultracompact Hii (UC Hii) regions (e.g. ISzymczak et al. 200^ 
and references therein). While methanol masers appeared to be 
associated with Hii regions at low angular resolution, subse- 
! quent high-resolution observations demonstrated that the masers 
' are generally not coincident with strong radio sources (e.g. 
Imishetal. 1998; Phillips et al. 1998). Further, statistics from 
blind surveys showed that most methanol masers did not have 
associated radio continuum, at least not at the few mJy level at 
a wavelength of 2 cm - only 10% of the sources detected in 
the Arecibo Methanol Maser Galactic Plane Survey (AMGPS; 
iPandian etal . 2007) had a counterpart at 21 cm within 23'' of 
the maser position, whi le two-thirds of the m ethanol masers in 
the Torun blind survey (ISzymczak et al.|[2Q02i) did not have a 6 
cm counterpart within V. 

The lack of radio continuum towards methanol masers has 
been interpreted in two ways - (i) The masers are associated 
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with massive star s prior to the origin of an ultracompact re- 
gion (IWalshet al.lil998.) . Rapid accretion onto a rnassive star 
can quench the formation of a H ii region (IWalmsleyll 19951) . (ii) 
Alternately, the exciting star has a mass lower than that of a pro- 
totypical high-mass star (spectral type B3 o r earlier), an d is thus 
too cool to produce a detectable H ii region (iPhillips et al.ll 19981) . 
iPhillips et al.l d 19981) concluded the latter based on the expected 
number of UC H ii regions in the Galaxy assuming that methanol 
masers overlapped with UC H ii regions for only 20% of their 
lifetimes. However, su bsequent work has favor ed the first hy- 
pothesis. For example, van der Walt et al.l (l2003l) argue that it is 
unlikely for the masers to be associated with stars of spectral 
type B5 based on the expected lifetime of maser emission and 
the number of methanol masers in the Galaxy. They further ar- 
gue that the eflficiency of exciting a maser would decrease dra- 
matically with decreasing stellar mass. 

Surveys for 6.7 GHz methanol masers towards sites of 
low- mass star formation have no t yielded any detections 
(Mi nieretal.1 120031: IXu et al.1 [2008h . Based on their survey, 
Mini er et aP (12003 ^ give a lower limit to the stellar mass as- 
sociated with bright (brightness temperature > 3 x 10^ K) 6.7 
GHz methanol masers as 3 M©. The bolometric luminosities of 
the sources associated with 6.7 GHz maser emission, estimated 
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using mid-infrared and IRAS data, were sugge st ive of stars 
with s trong ionizing radiation (IWalsh et al.lll999l) . IWalsh et aP 
(|2QQ1|) modeled the spectral energy distributions of the masers 
using IRAS and near-infrared data and found that sources for 
which spectral types of the central stars could be unambigu- 
ously identifi ed were consistent with massive, ionizing stars. 
iMinier et aP (2005) carried out a multi- wavelength study of five 
radio-quiet methanol masers, and found the environment around 
the masers to be characteristic of massive star forming clumps in 
earlier evolutionary phases than H ii regions. The observation of 
warm dust around metha nol masers with b olometric luminosi- 
ties greater than 10^ L© (Wals h et"aDl2003l) also suggested that 
the masers are associated with massive star forming regions. 

One of the limitations of these studies has been that the 
methanol maser sample is not complete. They are typically se- 
lected from a targe ted survey , and are hence potentially bi- 
ased. The work of lEUingsenl (l2006h . which investigated the 
mid-infrared colors of the masers using the Galactic Legacy 
Infrared Mid-Plane S urvey Extraordinaire survey (GLIMPSE; 
iBenja min et al.ll2QQ3l) , did use a complete sample, but no mod- 
eling of the spectral energy distributions (SEDs) was carried out 
on account of poly cyclic aromatic hydrocarbon (PAH) emission 
and silicate absorption lines falling in the Spitzer IRAC bands. 

Recently, a grid of 200,000 young s t ellar o bject (YSO) mod- 
els was developed by iRobitaille et al.l (|200o) spanning a wide 
range of evolutionary stages for different stellar masses to model 
the SED from optical to millimeter wavelengths. In addition, 
24 fim data is now available for the Galactic plane from the 
Spitzer legacy project, MIPSGAL (ICarev et al. 1 120091) . The 6'' 
beam of Spitzer at this wavelength and its much better sen- 
sitivity is a significant improveme nt to that of the Midcourse 
Source Experiment satellite (^MSX: llPrice et al.lll999l) . In this pa- 
per, we present multi-wavelength observations from centimeter 
to submillimeter wavelengths towards a sample of 20 6.7 GHz 
methanol masers selected from a blind survey. We combine our 
observations with existing data to determine SEDs from cen- 
timeter to near-infrared wavelengths. The combination of the 
radio continuum data and the results from SED modeling give 
very strong evidence that the methanol masers are associated 
with rapidly accreting massive stars. 



2. Sample, observations and data reduction 

2.1. The sample 

To determine the spectral energy distributions of 6.7 GHz 
methanol masers, we observed a sample of 20 masers with the 
Very Large Arra}Q (VLA) at 3.6 cm, L3 cm and 6.9 mm, the 
IRAM 30-meter telescope at L2 mm and the Ata cama Path 
Finder Experiment telescope, APE^ Glisten et al. 2006), at 870 
Jim. The source sample was selected from the AMGPS and in- 
cluded all sources detected between Galactic longitudes of 38.6° 
and 43.1°, and Galactic latitudes \b\ < 0.42°, the latitude con- 
straint being imposed for the sake of completeness. The peak 
flux density of the mas ers range from 0.2 Jy to 26 Jy, and their lu- 
minosities (taken from lPandian et al.,2009.) range from 9.9x10"^ 
to 2.6 X 10-5 Lq. 



^ The National Radio Astronomy Observatory is a facility of the 
National Science Foundation operated under cooperative agreement by 
Associated Universities, Inc. 

^ This publication is based on data acquired with the Atacama 
Pathfinder Experiment (APEX). APEX is a collaboration between 
the Max-Planck-Institut fiir Radioastronomie, the European Southern 
Observatory, and the Onsala Space Observatory. 



14 out of 20 sources have accurate positions (astrometry bet- 
ter than O.f ) determined from phase-referenced observations 
with the MERLIN interferometer (these observations were part 
of a larger project to determine accurate positions to the AMGPS 
sources and will be described in detail in Pandian et al. 2010, 
in preparation). The remaining sources have positions deter- 
mined from 24 yum point sources in MIPSGAL. It is generally 
found that 6.7 GHz methanol masers coincid e with 24 jim point 
sources with good accuracy (IXu etal .1120091 Pandian et al., in 
preparation). For instance, the mean deviation between accu- 
rate maser positions and their corresponding 24 yum MIPSGAL 
point sources for the 14 sources above is l.V\ Consequently, 
for sources without good interferometer positions, we can use 
MIPSGAL point sources to determine positions accurate to ~ 

r^ 

As indicated in iPandian & GoldsmithI (l2007h . the AMGPS 
positions have a root mean square positional accuracy to 7'\ 
which translates to an uncertainty of 18'' at the 95% confidence 
level. Hence, we looked for MIPSGAL point sources within 18" 
of the 6 methanol masers without good interferometer positions. 
In most cases, there was a single point source within the search 
radius that was defined as the counterpart for the maser. In cases 
where there were more than one 24 jim source, we looked at 
mid-infrared colors using the GLIMPSE point source catalog, 
and selected counterparts that had colors consistent with those 
of young stellar objects. 

The positions, peak flux densities, distances and luminosities 
of the masers in our sample are shown in Table 1 . 



2.2. Observations and data reduction 

The 3.6 cm observations were carried out on December 8, 2007 
with the VLA in the B -configuration. The correlator was set 
up in the standard continuum mode, and provided a total band- 
width of 100 MHz per polarization. The integration time was 2.5 
minutes per source. The source 3C286 was used as a primary 
and flux calibrator while J1856-h061 served as the phase cali- 
brator. The data were reduced using standard procedures in the 
Astronomical Image Processing System (AIPS) of NRAO. The 
full width at half-maximum (FWHM) of the synthesized beam 
was 0.9'' X 0.8", and the Icr noise in the images range from 0.13 
to 0.24 mJy/beam. 

The 1.3 cm observations were carried out with the VLA 
on March 22, 2008 in the C-configuration. The source 3C286 
was used as a primary flux calibrator while J1851-h005 served 
as the phase calibrator. The observations were carried out us- 
ing fast switching for good phase stability. Each fast switching 
cycle comprised of 140 s on the target followed by 40 s on the 
phase calibrator. The total on-source integration time was 6 to 13 
minutes per source. The data were reduced using AIPS. The Icr 
noise in the maps was around 0.12 mJy/beam, the synthesized 
beam being 1.0" X 0.9". 

A subset of the sources in the sample (see Table O were 
observed with the VLA at 6.9 mm in the D-configuration on 
September 2, 2008. As with the 1.3 cm observations, 3C286 was 
used as a primary flux calibrator while J1851-h005 served as the 
phase calibrator. Fast switching was used, each cycle comprising 
of 240 s on the target followed by 60 s on the phase calibrator, 
with the total on-source integration times ranging from 4 to 31 
minutes depending on the strength of the source. The data were 
reduced using AIPS. The root mean square noise level in the 
maps ranged from 0.4 to 1.7 mJy/beam and the FWHM synthe- 
sized beam sizes ranged from 2.7" x 2.2" to 3.8" x 2.3". 
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Table 1. 6.7 GHz methanol masers observed at multiple wavelengths to determine SEDs. The columns show the source name, equa- 
torial J2000 coordinates, peak flux density, distance and luminosity. Positions that are accurate to O.T' are derived from MERLIN 
observations, while positions accurate to V are derived from MIPSGAL 24 jim point sources. 



Source 



a 
(J2000) 



6 
(J2000) 



(Jy) 



d 
(kpc) 






38.66+0.08 


19 01 35.24 


05 07 47.4 


0.71 


16.3 + 0.9 


38.92-0.36 


19 03 38.66 


05 09 42.5 


1.29 


10.5 + 0.4 


39.39-0.14 


19 03 45.32 


05 40 42.6 


0.77 


4.3 + 0.5 


39.54-0.38 


19 04 52.6 


05 42 08 


0.24 


9.0 + 0.5 


40.28-0.22 


19 05 41.22 


06 26 12.7 


68.03 


4 Q+0.9 
•^ 6 


40.62-0.14 


19 06 01.59 


06 46 36.1 


6.97 


10.5 + 0.4 


40.94-0.04 


19 06 15.38 


07 05 54.5 


1.50 


10.0 + 0.5 


41.08-0.13^ 


19 06 48.5 


07 1100 


0.33 


8.4 + 0.6 


41.12-0.11 


19 06 50.2 


07 14 02 


0.63 


10.0 + 0.5 


41.12-0.22 


19 07 14.86 


07 1100.7 


1.23 


8.7 + 0.6 


41.16-0.20 


19 07 14.3 


07 13 19 


0.20 


8.7 + 0.6 


41.23-0.20 


19 07 21.38 


07 17 08.2 


7.55 


8.7 + 0.5 


41.27+0.37 


19 05 23.6 


07 35 06 


0.16 


11.5 + 0.5 


41.34-0.14 


19 07 21.84 


07 25 17.3 


25.40 


11.6 + 0.5 


41.58+0.04 


19 07 09.1 


07 42 25 


0.22 


11.5 + 0.5 


42.03+0.19 


19 07 28.18 


08 10 53.5 


30.10 


11.1+0.5 


42.30-0.30 


19 09 43.59 


08 1141.4 


5.07 


10.5 + 0.5 


42.43-0.26 


19 09 49.86 


08 19 45.4 


1.61 


7.9 + 0.8 


42.70-0.15 


19 09 55.07 


08 36 53.5 


4.03 


15.9 + 0.8 


43.08-0.08 


19 10 22.05 


08 58 51.5 


4.11 


11.2 + 0.5 



(1.3 + 0.2) X 10-^ 
(9.9 + 0.8) X 10-^ 
(9.9 + 2.3) X 10-^ 
(1.3 + 0.2) X 10-^ 

i.rjj-^xio-5 

(5.3 + 0.4) X 10-6 
(1.0 + 0.1) X 10-6 
(1.6 + 0.3) X 10-^ 
(4.4 + 0.5) X 10-^ 
(6.5 + 0.9) X 10-^ 
(1.0 + 0.2) X 10-^ 
(4.0 + 0.5) X 10-6 
(1.5 + 0.2) X 10-^ 
(2.4 + 0.2) X 10-^ 
(2.0 + 0.2) X 10-^ 
(2.6 + 0.3) X 10-5 
(3.9 + 0.4) X 10-6 
(7.0 + 1.4) X 10-^ 
(7.1+ 0.7) X 10-6 
(3.6 + 0.3) X 10-6 



Notes. (^^ An alternative position for this source based on MIPSGAL is (19^^ 06'" 49^0, T W 07'^. 



The 1.2 mm observations were carried out between January 
and March 2008 as part of "pooled obser vations" using the 1 17 
element MAMBO-II bolometer camera (iKrevsaet al.(ll998h at 
the IRAM 30-meter telescope. Each source was mapped using 
the standard on-the-fly technique such that a ~ 2' region was uni- 
formly sampled (although the actual extents of the maps are con- 
siderably larger). The time spent per map was around 12 minutes 
with pointing, focus and calibration scans being done at regular 
intervals. The data were reduced using the "mapCSF" procedure 
in the MOPSIC, which is part of the GILDAS software package 
(http://www.iram.fr/IRAMFR/GILDAS). The root mean square 
level in the maps was typically 8 to 10 mJy/beam with the 
FWHM beam size being ~ 10.5''. The maps were loaded into 
AIPS, and flux densities of point sources were extracted using 
the task "JMFIT". 

The 870 fim observations were carried out in December 
2007 using the La rge APEX Bolometer Camera (LABOCA; 
ISiringoet all 120091) at the APEX telescope. The field of view 
of the camera (~ IT) was mapped around each source (except 
for sources 41.08-0.13 and 41. 12-0.11, and 41. 12-0.22, 41.16- 
0.20 and 41.23-0.20, which are close enough to be covered in 
two separate fields of view) using raster scans. The time spent 
per map ranged from 5 to 22 minutes depending on the strength 
of the source and the weather conditions. The data were reduced 
using the Bolometer array data Analysis (Bo A; Schuller et al, in 
prep.) softw are. The data reductio n strategy was similar to that 
described by Schuller et al. (^20091) and involved flagging of bad 
pixels, removal of correlated noise, despiking, low-frequency fil- 
tering and first-order baselining of the timestream followed by 
building a map using appropriate weighting. These steps were 
carried out in several iterations, with a source model from the 
previous iteration being subtracted from the data, enabling more 
aggressive application of procedures such as despiking. The final 
Icr noise level in the maps ranged from 30 to 75 mJy/beam, the 
FWHM beam size being ~ 18.2''. As with the MAMBO data. 



flux densities of point sources were extracted using the AIPS 
task "JMFIT". 

In addition to our observations above, we also used existing 
data at mid-infrared wavelengths to determine the SEDs. 
While the point source catalogs of the GLIMPSE and 2MASS 
surveys provided data from 8.0 to 1.24 yum, the point source 
catalog of 24 yum sources from the MIPSGAL survey has not 
been released. Hence, we determined 24 yum flux densities 
directly from the mosaics released by the MIPSGAL team. 
To measure the flux densities, we followed the procedure 
outlined in the cookbook in the Spitzer Science Center website 
(http://ssc.spitzer.caltech.edu/dataanalysistools/cookbook/25/). 
For sources that were saturated in the MIPS images, we used 
flux densities fro m the Galactic Pl ane survey carried out with 
the MSX satellite ( JEgan et al.ll20Q3h . 



3. Results 

No continuum emission was seen in 14 out of 20 sources at both 
3.6 cm and 1.3 cm. At the native resolution, six sources were 
detected at 1.3 cm, of which five were detected at 6.9 mm and 
two were detected a t 3.6 cm and at 6 cm in the CORNISH survey 
d Purcell et al.ll2008h . The latter two sources have spectral indices 
consistent with optically thin free-free emission, while the three 
sources detected at 6.9 mm have spectral indices greater than 
0.3, which is consistent with optically thick free-free emission. 
The nature of the source 42.70-0.15, that was detected only at 
1.3 cm is not clear. However, its 3cr flux density limit at 3.6 cm 
would be consistent with optically thick emission between these 
wavelengths. 

One of the concerns in the above analysis is the presence of 
extended emission that is resolved in our observations. Several 
UC H II regions have been observed to b e embedded in extended 
ra dio continuum eniissio n (JKurtz et al.| [T999). The recent work 
of lLongmore et al.l (l20Q9h also found a significant fraction of 6.7 
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GHz methanol mas er sources that had no 3.6 cm radio contin- 
uum as measured by lWalsh et al.l (Il998h . to have extended emis- 
sion with the total flux at times being several orders of magni- 
tude discrepant from previous limits. Missing an extended radio 
continuum component due to spatial filtering of an interferome- 
ter can lead to errors in the classification and inferred age of the 
source. 

We hence re-imaged all sources at 1.3 cm by tapering the 
uv-coverage of the data. This decreases the angular resolution 
of the data (at the expense of slightly increased noise in the im- 
ages) enabling identification of extended emission. The 1.3 cm 
observations are sensitive to a largest angular scale of ~ 3O'0. 
With the exception of 39.39-0.14, no extended emission was 
found in the sources. Although this does not rule out the pres- 
ence of extended emission on scales larger than 30'', it sug- 
gests a lack of continuu m sources similar to those discovered by 
iLongmore eTaP (12009 ). In the case of 39.39-0.14, the extended 
emission led to the discovery of a counterpart at 3.6 cm leading 
to re-classification of the nature of radio continuum emission 
from optically thick to optically thin emission between 3.6 cm 
and 1.3 cm. A more detailed discussion of this source can be 
found in gni 

At 1 .2 mm wavelength, 1 8 out of 20 sources were detected as 
point sources, of which 17 were also detected at 870 yum. The 1.2 
mm maps of the sources are shown online (Fig.|4]). One source, 
41.16-0.20, was found to be in the middle of a filament at both 
1.2 mm and 870 /im, without any point source to extract flux 
densities from. All sources had counterparts at 24 jjm (see ^2.11) , 
and in the GLIMPSE point source catalog with flux densities in 
at least two IRAC bands. In addition, 7 sources also had coun- 
terparts in the 2MASS point source catalog. The flux densities 
of the sources from centimeter to submillimeter wavelengths are 
shown in Table [2l while the infrared flux densities are shown in 
Tables 



3.1. Notes on selected sources 

38. 66+0. 08 - The 3.6 and 1 .3 cm images show a prominent point 
source along with some hint of extended emission. The point 
source flux density measured with "JMEIT" is 14.6, 11.7 and 
11.9 mJy respectively at 3.6, 1.3 and 0.69 cm suggestive of a 
flat spectral index. However, the 6 cm flux density measured by 
CORNISH is only 6.97 mJy, and the resulting spectral index (de- 
fined as 5'y oc y^ where 5'y is the flux density at frequency v and 
a is the spectral index) between 6 cm and 3.6 cm is -hi. 3 which 
is consistent with optically thick free-free emission. Using the 
Altenhofl^ approximat ion, the free-free optical depth is given by 
(lAltenhoffetal.1119601) 



Tv = 



0.082 



/"• 



dl 



(1) 



where v is in GHz and J n^ dl, the emission measure, is in pc 
cm"^. Since the free-free emission in this source is optically 
thick at 6 cm (y = 4.86 GHz), but optically thin at 3.6 cm (y 
= 8.46 GHz), the emission measure is between 8.5 x 10^ and 
2.7 X 10^ pc cm"^, and can be classified as a UC H ii region. 

It is to be noted that the 1.4 GHz NRAO VLA Sky Survey 
(NVSS: [Condon et al.1 [19981) reveals a 165 mJy source about 
8.5'' away from the maser position, with a spatial extent of 
63'' X 58". Although this is well above the maximum angular 




http://www.vla.nrao.edu/astro/guides/vlas/current/node 1 0.html 



19^03^46" 44S 42S 4QS ggS 3gS 34S 

Right Ascension (J2000) 

Fig. 1. The 870 jim LABOCA image of source 38.92-0.36 is 
shown in color, overlaid on which is the 1 .2 mm MAMBO image 
in green contours. The contours start at 10% of the peak emission 
(550 mJy/beam) in steps of 5% thereafter. The black cross shows 
the location of the methanol maser. 



scale that can be recovered by our observations, we do not re- 
cover flux above the point source component (which is coinci- 
dent with the maser position) indicated in Table [2] when imag- 
ing the source with a tapered uv-coverage. Considering that this 
source is at a distance of 16.3 kpc from the Sun, the NVSS re- 
sults are indicative of a star forming complex with a classical H 11 
region, embedded within which is the UC H 11 region detected in 
our observations. 

The MAMBO observations show a complex structure with 
significant extended emission (which is consistent with the ex- 
tent of the NVSS emission) that cannot be fully recovered from 
the small size of the map, and consequently no 1 .2 mm flux den- 
sity is reported in Tabled However, the LABOCA map does not 
show any extended emission within the map sensitivity. 

38.92-0.36 - The LABOCA data shows a strong point 
source at J2000 coordinates (19^ 38^^ 38^9, 5° 10' 0") with a 
bow- shock like morphology (Fig. [T]). No point source is seen 
at the location of the maser, and consequently, no 870 yum flux 
density is reported in Tabled The higher resolution MAMBO 
data resolves the single LABOCA source into multiple sources 
including a point source at the maser site. 

39.39-0. 14 - At native resolution, this source appears to be a 
hypercompact H 11 (HC H 11) region since there is no counterpart 
detected at 3.6 cm, while the spectral index between 1.3 cm and 
6.9 mm suggests optically thick free-free emission. However, 
imaging the source with a tapered uv distribution shows the pres- 
ence of extended emission with a total 1.3 cm flux of 3.0 mJy, 
the extended component accounting for 50% of the flux. A 3.6 
cm image made with a similar tapered uv-coverage shows a 3.5(T 
point source with a total flux of 3.5 mJy, indicating that the emis- 
sion is optically thin between 3.6 cm and 1.3 cm. There is no 
de tection of the source in the CORNISH survey or the catalog 
of lBecker et al.l (Il99ll) which may be due to a combination of 
sensitivity and spatial filtering by the interferometer. Since the 
emission measure is less than 2.7 x 10^ pc cm"^, this source is 
classified as a UC H 11 region. The uv-tapered images also led 
to the identification of a nearby UC H 11 region (J2000 coordi- 
nates: 19^ 03'^ 46^0, 5° 40' 42") with 3.6 cm, 1.3 cm and 6.9 
mm fluxes of 4.3, 3.3 and 3.4 mJy respectively. 

40.28-0.22 - The 1.3 cm flux density in this source is 
0.4 ± 0.2 mJy, while the 6.9 mm flux density is 6.0 ±1.2 mJy, 
which gives a spectral index of 4. 1^^^, which is inconsistent with 
free-free emission. This suggests that dust emission contributes 



J. D. Pandian et al.: Spectral Energy Distributions of 6.7 GHz methanol masers 

Table 2. Centimeter to submillimeter wave flux densities measured towards the 20 6.7 GHz methanol masers. The indicated upper 
limits are 3cr and all flux densities are in mJy. 



Source 


5(3.6 cm) 


5(1.3 cm) 


5 (6.9 mm) 


5 (1.2 mm) 


SiSlOjim) 


38.66+0.08 


14.6 + 0.3 


11.7 + 0.2 


11.9 + 3.6 


a 


1390 + 320 


38.92-0.36 


<0.6 


<0.3 




1210+130 


b 


39.39-0.14 


3.5 + 1.4 


3.0 + 0.5 


5.0+1.3 


1514 + 48 


4550 + 170 


39.54-0.38 


<0.7 


<0.3 


< 1.8 


154 + 62 


480 + 120 


40.28-0.22 


<0.8 


0.4 + 0.2 


6.0+1.2 


3156 + 35 


12520 + 160 


40.62-0.14 


<0.8 


0.9 + 0.3 


3.9 + 1.7 


1657 + 56 


6280 + 280 


40.94-0.04 


<0.7 


<0.3 




53 + 26 


170+110 


41.08-0.13 


<0.6 


<0.4 




680 + 120 


2670 + 300 


41.12-0.11 


<0.5 


<0.4 




390 + 52 


1360 + 230 


41.12-0.22 


<0.6 


<0.4 




367 + 55 


3130 + 310 


41.16-0.20 


<0.5 


<0.4 




c 


c 


41.23-0.20 


<0.5 


<0.4 




264 + 35 


1590 + 250 


41.27+0.37 


<0.5 


<0.4 




309 + 60 


860 + 100 


41.34-0.14 


<0.5 


<0.4 




327 + 39 


1140 + 230 


41.58+0.04 


<0.4 


<0.4 


< 1.3 


<28 


< 110 


42.03+0.19 


<0.4 


<0.4 




572 + 81 


1590 + 260 


42.30-0.30 


<0.4 


<0.4 


< 1.5 


538 + 57 


2610 + 200 


42.43-0.26 


59.5 + 1.2 


49.4 + 0.8 


45.9 + 5.5 


670 + 36 


2230 + 160 


42.70-0.15 


<0.4 


0.8 + 0.2 


<4.5 


177 + 44 


500 + 130 


43.08-0.08 


< 1.5 


<0.5 




595 + 67 


2220 + 49 



Notes. ^^^ 1.2 mm flux density uncertain due to presence of extended emission that is not fully recovered. 
^^^ Source not detected at 870yL/m due to poorer resolution of LAB OCA. 
^'^^ No point source seen due to source being in the middle of a filament. 

Table 3. Infrared flux densities of the 20 6.7 GHz methanol masers. 



Source 


MIPS 




MSX 






IRAC 






2MASS 






24yt/m 


21.3 //m 


14.7 yum 


12.1 yum 


8.3 yum 


8yum 


5.8 yum 


4.5 yum 


3.6 yum 


K 


H 


J 




(Jy) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


38.66+0.08 


2.37 


2.37 






0.14 




19.5 


7.70 


0.80 








38.92-0.36 


0.051 














3.23 


0.40 








39.39-0.14 


3.88 










113 


144 


72.1 










39.54-0.38 


0.13 










4.74 


4.89 


4.35 


1.44 








40.28-0.22 


>5.9^ 


7.15 


1.03 


0.064 




289 


225 


252 


129 


19.2 


3.70 




40.62-0.14 


a 


29.1 


4.94 


2.29 


0.98 


1522 


962 


355 




3.03 






40.94-0.04 


0.085 














0.87 


0.42 








41.08-0.13A 


0.055 










6.12 


8.11 


7.65 


2.51 








41.08-0.13B 


0.081 














4.22 


2.38 


2.52 


1.56 




41.12-0.11 


4.03 


2.73 


2.73 


1.83 


1.27 


1482 


1526 




317 


6.02 






41.12-0.22A 


0.15 














3.87 


2.02 








41.12-0.22B 


0.15 












1.23 


1.34 


0.32 








41.16-0.20 


0.17 












1.23 


1.34 


0.32 








41.23-0.20 


0.26 


1.13^ 


0.97^ 


1.06^ 


0.69^ 






8.07 


4.37 








41.27+0.37 


2.33 


1.51 


0.81 


0.95 


0.37 


343 


183 


106 


54.7 


7.88 






41.34-0.14 


0.91 










20.5 


11.9 


5.59 










41.58+0.04 


0.019 














0.55^ 


0.73^ 








42.03+0.19 


>5.3^ 


15.96 


9.40 


5.97 


2.42 


1796 


2362 


681 


463 


38.8 


2.84 




42.30-0.30 


2.49 










89.6 


47.9 


27.2 


6.58 








42.43-0.26 


a 


67.58 


18.84 


12.65 


4.47 




1118 


991 


565 


210 


60.0 


8.01 


42.70-0.15 


2.88 










36.7 


11.9 


3.38 


1.27 








43.08-0.08 


0.83 












3.24 


5.02 


1.06 









Notes. ^^^ Saturated in MIPS data. 

^^^ Significant uncertainty in the MSX flux densities on account of the extended structure seen in the higher resolution MIPS data. 

^^^ An offset of 2.8'' between the GLIMPSE and MIPSGAL point sources casts some doubt about the association of the two sources. 



to the 6.9 mm flux density. One can obtain an estimate of the con- 
tribution of dust emission by fitting the 1.2 mm, 870 yum and 24 
yum flux densities with a grey-body of the form By(Td) (1 - e''^") 
where T^ is the dust temperature and Zy is the optical depth and is 
given by to(v/vo)^, where yS is the dust emissivity a nd is assumed 
to be 2. It should be noted that previous work (e.g. lMinier et al.l 



12005^ typically requires more than one temperature component 
to fit the SED at wavelengths shorter than ~ 50 yum. However, 
since we have only three data points, it is not possible to carry 
out a multi-temperature fit in our case. The consequence of us- 
ing a single temperature model would be to overestimate the dust 
temperature and consequently the contribution of dust emission 
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Right Ascension (J2000) 






Fig. 2. The 1.3 cm VLA image of the source 40.62-0.14. The 
red cross marks the location of the methanol maser. 



to the 6.9 mm emission. We obtain a dust temperature of 47 K 
from this fit, and estimate the contribution to the 6.9 mm flux 
density to be 3.4 mJy. Taking this into account, the free-free 
spectral index is 2.8 ± 1.6, which is marginally consistent with 
optically thick free-free emission. 

40.62-0.14 - The 1.3 cm map of this source (Fig. [2]) shows 
two sources, one of which is coincident with the maser location. 
The continuum associated with the maser source is not detected 
at 3.6 cm, and has 1.3 cm and 6.9 mm flux densities of 0.9 ± 0.3 
mJy and 3.9 + 1.7 mJy respectively. The other stronger source, 
at J2000 coordinates (19^ 06^^ OP .48, 6° 46' 35''.5), has 3.6 cm, 
1.3 cm and 6.9 mm flux densities of 3.0 ± 0.6 mJy, 2.0 ± 0.2 mJy 
and 3.8 ± 1.3 mJy respectively. The two sources are unresolved 
in the MAMBO and LABOCA maps. A grey-body fit to the 1.2 
mm, 870 yum and 21 jim data (the latter being from MSX since 
the MIPS data is completely saturated) gives the contribution of 
dust to the 6.9 mm emission as 1.7 mJy. Taking into account the 
error bars, the maser source is consistent with a HC H ii region, 
while the nearby source appears to be a UC H ii region. The latter 
is not detected in CORNISH (3cr limit of 0.9 mJy), which sug- 
gests that the emission measure is between 8.5 x 10^ and 2.7 x 10^ 
pc cm-^ similar to that of 38.66+0.08. The GLIMPSE source in 
Table [3] is coincident with the maser source. 

41.08-0. 13- There are two possible MIPSGAL counterparts 
for this source as indicated in Table [T] The infrared flux densities 
of the two sources are given in Table [3] by the suffixes 'A and 
'B' respectively ('A referring to the position in Tabled and 'B' 
referring to the source in the footnote). Based on the properties 
of the GLIMPSE source, the source 'A is more likely to be the 
infrared counterpart for the maser. 

41.12-0.22 - There are two possible GLIMPSE point 
sources (labeled 'A and 'B' in Table [3]) that could be mid- 
infrared counterparts to the maser. The separation of the 'A and 
'B' sources from the maser are 0.7'' and 1.4" respectively, while 
the separation between the two sources and the 24 yum MIPS 
source are 2.0" and 2.1" respectively. 

42.43-0.26 -This is the strongest source at centimeter wave- 
lengths in our sample, with a 3.6 cm flux density of 59.5 mJy. 
The 6 cm flux density measured by the CORNISH survey is 
105.4 mJy, which is much higher than what is expected from ex- 
trapolating the 3.6 cm flux density using optically thin free-free 
emission. However, this is most probably due to the complex and 
extended structure in the source which can be seen in the 3.6 cm 
and 1.3 cm maps. Since the emission is optically thin to at least 



6 cm, the emission measure is lower than 8.5 x 10^ pc cm"^, and 
is a UC H II region. The relatively evolved nature of the source is 
also evident in this being the only source detected in the 2MASS 
survey at /-band. 

42. 70-0. 15- This source is detected at 1 .3 cm at the 8(T level 
though the signal to noise ratio in the integrated flux density is 
less. The 3.6 cm image shows a 3cr source with a peak inten- 
sity of 0.4 mJy/beam located 0.6" away from the 1.3 cm source. 
However, it is not clear whether this source is real, and it is not 
possible to obtain a good estimate of the integrated intensity of 
the source. Consequently, this is reported as a non-detection at 
3.6 cm in Table O Assuming that the 3.6 cm source is real and 
is unresolved, the source would be classified as a hypercompact 
H II region. 

4. Discussion 

Using eq. (1) for the optical depth of free-free emission under 
the Altenhoff' approximation, an optical depth of unity at 1.3 cm 
and 6.9 mm correspond to emission measures of 2.1 x 10^ and 
8.4 X 10^ pc cm"^ respectively. Hence, as mentioned in ^3.11 
the two sources that have continuum emission that is optically 
thick between 1.3 cm and 6.9 mm (40.28-0.22 and 40.62-0.14) 
can be classified as hypercompact Hii regions. Three sources, 
38.66+0.08, 39.39-0.14 and 42.43-0.26, are ultracompact Hn 
regions, while the source 42.70-0.15 is potentially a hypercom- 
pact H II region. 

The lack of centimeter wave counterparts for the major- 
ity of the sample is consistent with the results from previous 
work. This could arise from two scenarios: the methanol masers 
could be associated with very early phases of massive star for- 
mation that had their H ii regions quenched by rapid accretion. 
Alternatively, the masers could arise around intermediate mass 
stars whose ionizing flux is too weak to result in detectable H ii 
regions. 

To try to distinguish between the two scenarios, we adopt 
two approaches. First, we use the centimeter wave data to obtain 
upper limits on the sizes of any undetected H ii regions and the 
ionizing radiation of the central star. Assuming that any unde- 
tected H II region is optically thick at frequency, v, and that the 
electron temperature in the Hii region is Tg = 10"^ K, the solid 
angle of the source, fig is related to the flux density, Syhy 



Sy - By(Te)^s 



(2) 



where By(T) is the black body function. Using the 3cr flux den- 
sity limits at 1.3 cm (y = 22.46 GHz), we calculate the limits on 
the source size, which are tabulated in Table IH Further, we have 
a constraint on the emission measure, which to first order can be 
taken to be ri^l where rie is the electron density and / is the size of 
the source. We then use ionization equilibrium to calculate limits 
on the ionizing flux of the central star, the results being tabulated 
in Table m 

The limits on the ionizing flux from the central object are 
not inconsistent with the presence of massive stars since the op- 
tical depth of free-free emission at 1.3 cm (tu in Tabled is 
unknown. However, the limits on the source size show that any 
ionized region must be confined to a volume that is comparable 
to that of the solar system (up to the Kuiper belt). As mentioned 
previously, these results could also be explained by the central 
objects being intermediate mass stars. 

It is worth pointing out that the HC Hii regions associ- 
ated with CRL 2136, ^ V33 A , NGC 2591, NGC 7538 IRS 9 
(iMenten & van der TakI 120041: Ivan der Tak & MentenI l2005h . 
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Table 4. Constraints on the sizes of ionized regions and ionizing 
flux from the central star for sources with no detections at 1.3 
cm. The columns show the source name, the limits on the source 
size, / in milliparsecs and AU, and the limit on the ionizing flux 
divided by the optical depth of free-free emission at 1.3 cm. 



Source 


/ 


/ 


M/Ti.3 




(mpc) 


(AU) 


10^5 photons s-i 


38.92-0.36 


0.47 


98 


2.5 


39.54-0.38 


0.41 


84 


1.9 


40.94-0.04 


0.45 


93 


2.3 


41.08-0.13 


0.38 


78 


1.6 


41.12-0.11 


0.45 


93 


2.3 


41.12-0.22 


0.39 


81 


1.7 


41.16-0.20 


0.39 


81 


1.7 


41.23-0.20 


0.39 


81 


1.7 


41.27+0.37 


0.52 


107 


3.0 


41.34-0.14 


0.52 


108 


3.1 


41.58+0.04 


0.52 


107 


3.0 


42.03+0.19 


0.50 


103 


2.8 


42.30-0.30 


0.47 


98 


2.5 


43.08-0.08 


0.51 


104 


2.9 



which all have bona fide associated class II CH3OH masers, 
would not or only marginally have been detected at the flux den- 
sity limits of our VLA observations. Consequently, the sizes that 
these authors determine for the above sources (80, 150, 20, and 
50 AU, respectively) are all at or lower than the upper limits on 
sizes we present in Tabled 

An alternate approach to the problem is to fit the infrared 
flux densities using models of young stellar objects which pro- 
vides constraints on t he physical parameters of the sources. We 
used the SED fitter of iRobitaille et al.l (l2007l) , which uses a grid 
of 200,000 precomputed model SEDs spanning a wide range of 
evolutionary stages for diff'erent stellar masses (assuming that 
stars form by accretion through a disk and envelope). The fitting 
procedure involves interpolation of the model fluxes to the aper- 
tures used to perform the photometry, scaling them to a number 
of distances between Jmin and Jmax (calculated from the uncer- 
tainties in the distance determination), followed by fitting to the 
data with the visual extinction, Ay, being a free parameter (al- 
lowed to vary between and 100 magnitudes in our fitting). The 
extinction model used by Robi taille et al.l (l2007h is the method of 
iKim et al . (1994 ) modified for the niid-in frared extinction prop- 
erties derived by llndebetouw et aU (l2005h . 

To account for possible calibration errors between diff'erent 
data sets, we set the minimum uncertainty in the flux densities 
to be 10%. For the millimeter and submillimeter data, the aper- 
ture sizes were set to the source sizes measured by "JMFIT". The 
measured source sizes for the 1 .2 mm data typically ranged from 
10 to 40'', while for the 870 yum data, the sizes were 20 to 40''. 
The visual extinction was set to vary between and 100 mag- 
nitudes. The SED fitter provides the parameters of the models 
that fit the data in order of increasing ;^^ goodness of fit. We cal- 
culate the mean and standard deviation of the stellar mass (M^), 
temperature (T^), radius (R^), envelope accretion rate (M), age, 
total luminosity (Ltot), visual extinction (Ay), and inclination to 
the line of sight (/) from the models using weighting in accor- 
dance with the x^ probability distribution. The standard devia- 
tion is typically asymmetric about the mean, and in some cases 
unrealistically small when few models fit the data and the dif- 
ference in x^ between successive models is large. The results 
of this study are summarized in Table \5\ and the model fits are 
shown in Figure [3l We do not show any disk parameters since 



the models are dominated by the envelope and the disk mass is 
small compared to the mass of the central star. 

Three sources are not included in Table \5\ The sources 
41.16-0.20 and 41.58+0.04 do not have 1.2 mm or 870 yum flux 
densities, and in the absence of at least one of these two flux den- 
sities, SED fits are poorly constrained. In addition, a large num- 
ber of models fit the data for source 40.94-0.04, and so mean- 
ingful parameters cannot be deduced by this method. As noted 
in ^3. H and Table [3l there are two possible infrared counterparts 
for 41.08-0.13 and 41.12-0.22. Hence, both possible SEDs have 
been modeled and the results indicated by suffixes 'A and 'B' in 
Table [5] and Figure O 

It is to be noted that uncertainties shown in Table [5] only take 
into account the x^ statistics of multiple models fitting the data, 
and do not show the systematic uncertainties from the underly- 
ing model assumptions. For example, model parameters such as 
stellar mass and age are derived from the lu minosity and teni - 
perature of the star using evolutionary tracks (lRobitailldl2QQ8l) . 
Consequently, uncertainties in the evolutionary tracks translate 
to uncertainties in the masses and ages. These uncertainties have 
not been quantified and are not shown in Tabled Further, the en- 
velope accretion rate shown in Table [5] is derived from the den- 
sity profile measured from the SED through the assu mption of 
a free-fall rotational collapse model (lRobitailldl2008l) . Since the 
millimeter and submillimeter data (which primarily trace the en- 
velope) have relatively poor angular resolution, source multiplic- 
ity potentially leads to over-estimation of the envelope accretion 
rates. 

Keeping the caveats above in mind. Table \5\ shows that the 
YSO models that fit the infrared and submillimeter data points 
are rapidly accreting massive stars. The central stars have masses 
larger than ~ 8 M©, and considering the rapid accretion, the final 
mass of the central star is likely to be much higher. In the case of 
the source 41.58+0.04, the upper limit at 1.2 mm corresponds to 
an isothermal mass of between 38 and 87 M© depending on the 
prescription of du st opacity (e.g. lOssenkopf & Hennin gill 9941: 
iHildebrandl 1 1 9 8 3l) assuming the dust temperature to be 30 K, 
which is adequate to host a massive star (although the constraint 
is much tighter if one uses the Icr rather than 3cr limit). 

The accretion rates and stellar temperatures in Table \5\ give 
further confirmation of the results at centimeter wavelengths. 
The stellar temperatures are mostly too cool to form H 11 regions 
(which is consistent with their young age and rapid accretion; 
e.g. Hosokawa & Omuk ai 2009). Moreover, as pointed out by 
IWalmsleyl (Il995l) . the "critical accretion rate" required to con- 
fine the Hii region to a volume close to the stellar surface is 
10""^ to 4 X 10"^ Mq yr"^ for stars of spectral type 05 (~ 60 M©) 
to BO (~ 17 Mq) respectively. The accretion rates deduced from 
the SED fitting are well above the critical accretion rate. This 
result is unlikely to change even accounting for source multi- 
plicity, which is common in the context of massive star forming 
regions. However, this calculation assumes that the accretion is 
spherical, which is not true in practice. This might explain why 
some of the sources are seen as ultracompact and hypercompact 
H II regions in spite of having accretion rates that would prohibit 
such regions from forming. 

We cannot deduce a minimum mass for the central objects 
associated with 6.7 GHz methanol masers since our sample size 
is relatively small. Our work does not rule out the possibility 
of a small number of masers being associated with stars with 
spectral type later than B3. However, accepting the assump- 
tion that all the methanol masers in our sample are massive 
YSOs (the mystery of 41.58+0.04 not withstanding), it is in- 
teresting to see that only six sources out of twenty are associ- 
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Fig. 3. Fits to the SEDs of 6.7 GHz methanol masers. The solid line shows the best fitting model, while the grey lines show other 
models whose {x^ -;^^^^^)/datapoint < 3. The dashed line shows the stellar photosphere corresponding to the central star in the best 
fitting model in the absence of circumstellar dust, but including interstellar extinction. 
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Fig. 3. Continued. 
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Table 5. Physical parameters of the young stellar objects associated with 6.7 GHz methanol masers from modeling their SEDs. The 
columns show the source name, age, mass, temperature and radius of the central object, envelope accretion rate, total luminosity, 
visual extinction to the source and model inclination to the line of sight. Note that the uncertainties are calculated solely from the;^^ 
statistics of different models that fit the data and do not take into account uncertainties from source multiplicity or other systematic 
uncertainties in the underlying model assumptions. 



Source 



Age 
(10^ yr) 



Q C+23.3 
^•-^-0.3 

4 9+0 

^^-11 

27+2 

30+f 

-^•^-2.3 

6.5+1-6 

7 4+23.4 
' -^-2.8 

18+15 



(Mo) 



16.8!^ 

99 7+0!2 

10.9!^:^ 
^•^^ 

38.4+^2 
19.4+^-2 

16.4+3-8 
27 1 +0-2 

97 Q+11.1 

21.6+1-4 
18.8 + 3.3 

12.2+6-^ 
16 2+^-^ 
11 5+4-^ 
21.3+^-0 
13 2+3-2 
16.4+^-2 
15.7+2-0 
24 5+2-9 



(K) 

"8600™^ 

8900+100 

13 300+1200 
4600:^00 

15 400+0 

32 300:^300 

5300+0 
4200+^05 

14 000+2800 

07^0+13 100 
y I JV7_24QQ 

16000+8^00 

4400+300 
4300+J^^ 

49QQ+8400 

36 300"+?^oo 
4300+8800 

30900+2100 
10 900+4600 



(Ro) 



62+1 

87+f 

18!20 
38!i5 
65t 

y+O 

290!Oo 
100+^™ 

110+31 

140:J 

o^_44 

"-0 
120^84 

5.8+O2 

41+31 
140+140 

^^ -100 



M 

(10-3 Mo/yr) 



(103 Lo) 



6.6!^2 

y y+O.l 

2.6+ 

1.7+ 



42^10 

Q 1 +0.8 
51 + 1.72 

2io:?o 

17+10 
-.0 

23+3 

95!?9 
38!^ 

34!f 

4 7+22.3 

5 4+0.6 

4 0+13-4 

-a.v 

55+22 

29!g 

17+23 

6r38 



(mag) 



55^^ 
26^48 

21+4 
26+5 



31!? 



37+2 
17+6 

15+9 
50+10 

7 + 1 

34+9 

24+1 
-0 

26+2 

10!S 

66!!^ 

27^8 



n 



38.66+0.08 

38.92-0.36 

39.39-0.14 

39.54-0.38 

40.28-0.22 

40.62-0.14 

41.08-0.13A 

41.08-0.13B 

41.12-0.11 

41.12-0.22A 

41.12-0.22B 

41.23-0.20 

41.27+0.37 

41.34-0.14 

42.03+0.19 

42.30-0.30 

42.43-0.26 

42.70-0.15 

43.08-0.08 



14 

2.8:}-0 

9 7+2:0 
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ated with HC Hii and UC Hii regions (although the absence 
of HC Hii region s could be due to sensitivity limitations as 
discussed earlier). IWalsh et aH (11998) suggested that methanol 
masers are destroyed as a UC Hii re gion develops. How ever, 
there are sources such as W3(QH) (JMentenet al.l 119921) and 
NGC 7538 IRSl (' Minier et alJ I2QQ0I: iMoscadelli et al.1 l2QQ9h 
where methanol masers are associated with strong and well 
developed UC Hii regions (sizes > 1000 AU). Furthermore, 
iHartquist etaP ([1995.) argue that OH and CH3OH can coexist 
in the envelope of a UC H 11 region. 

One might argue that the OH maser switches on at a later 
evolutionary stage and lasts longer. This appears to observa- 
tionally confirmed: many methanol masers do not have associ- 
ated OH masers. Although this may, again, be a selection eff'ect 
caused by the 6.7 GHz methanol maser's intensity to be on av- 
erage ^ 7 times high er than the strongest transition of the OH 
maser (iCaswell et al.i.l99 5), it might indeed be the case. Both 
the high OH and CH3OH abundances are the result of frozen out 
OH and H2O being evaporated off' dust grain mantles by heating 
from the igniting pro to star. However, the production of OH by 
photodissociation of H2O requires ultraviolet radiation, which 
necessitates the presence of an UC H 11 region. 

In addition, the estimated lifetime of the m ethanol maser 
emission is only a few x 10"^ yr Ivan der W alt' 2005). These would 
argue against the suggestion that the development of the UC H 11 
region destroys the methanol maser emission. However, our data 
does show that the likelihood of the presence of maser emission 
decreases as the massive star evolves into an H 11 region. Thus, 
our work would be consistent with a statistical statement that 
6.7 GHz methanol masers are most likely to be associated with 
massive stars prior to the formation of a HC H 11 region, and dis- 
appear during the UC H 11 phase. 

As mentioned previously, one of the big caveats in this study 
is the sources may resolve into multiple components at high an- 
gular resolution. This is especially so for the millimeter and sub- 



millimeter data which have much poorer resolution compared 
to other data. While the fitting procedure takes into account the 
aperture used for calculating the ffux density of the source, the 
assumption is that the dusty envelope seen in the millimeter and 
submillimeter data harbors a single massive YSO. Data at high 
angular resolution is required to determine the severity of the 
eflfect of multiplicity in these results. Observations with exist- 
ing and upcoming submillimeter interferometers will be able to 
address this question. A more complete description of the lim- 
itations of the SEP modeling employed here can be found in 
lRobitailld(l2008i) . 



5. Conclusions 

We have obtained SEDs of a complete sample of 20 6.7 GHz 
methanol masers from centimeter to near-infrared wavelengths 
by combining new observations with existing data. Centimeter 
wave continuum is not seen in 14 sources, while three sources 
each are associated with HC Hii and UC Hii regions respec- 
tively. The upper limits on the sizes of any ionized region around 
sources with no radio continuum are consistent with such re- 
gions being confined close to the stellar surface. Modeling of the 
SED from millimeter to near-infrared wavelength strongly sug- 
gests that the young stellar objects associated with the masers are 
rapidly accreting young massive stars. Although the relatively 
small size of our sample does not allow us to put any constraints 
on the minimum mass of the central star of the YSOs associated 
with the masers, this work shows that most 6.7 GHz methanol 
masers are associated with massive stars prior to the formation 
of a hypercompact H 11 region. 
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Fig. 4. Continued. 



J. D. Pandian et al.: Spectral Energy Distributions of 6.7 GHz methanol masers, Online Material p 3 



4L58 + 0.04 



42.03+0.19 



+ 7°43'00" - 









I 






-50 



10^ 



08^ 



Right Ascension (J2000) 
42.30-0.30 




o +8°12'00' 
o 



o 



f36'30' 



19^09^58" 56' 54' 52' 

Right Ascension (J2000) 



29' 28' 27' 26' 

Right Ascension (J2000) 
42.43-0.26 




o +20'00' 





IT- 1 


_ 




^ 


^ra 


^" 




■ 


r 








, ■ 










19^09"^46' 44' 42' 

Right Ascension (J2000) 

42.70-0.15 


^^^■H ^ ^21 


■ 


■ 






■ ^M ^ 






I 


^^J 


r 


f" "i 




-1 


^ 1 


1 


i y 




■ 


f'- 


1 


Mr 


■ 


1 








^^1 


m H 


■ I 


\m 


LH 


1 


■^^B 




■ ■ 


^^^^^^^ 


■ 


^H 


1 


1 




mm m 




^^ 


■ - 


J 




L -^ 


■ 


■ 




H 


— ■ 


I 






H 




\ 


H 


a 


■ ^ 


— 1 


1 


■ 






■ 




^ 


■1 




I 






■ 


I 


wm 






1 


■ ■ 


r J 


^■~ 


1 


J 


■ 


^1 




■ 1 






J 


^ 


■ 


L 






L 






^^^H 




■1 ■ n 


■ 


!■ 




^ 

>. 



19^09"'52' 51' 50' 49' 

Right Ascension (J2000) 
43.08-0.08 



60 o 

o 

B ^ 

\ .2 
- -I20 B .B 



i 



P3H 

8°59'30"L TlV^^^^^^I 


^p 


L^T^ 


Ml 


+ 58'30"H^LJP! ■ 


^ 


CirtA 


M 



- - 100 






I 



19^10™26' 24' 22' 20' 

Right Ascension (J2000) 



Fig. 4. Continued. 



